Graphite in the bi-layer regime: in-plane transport 



D. B. Gutman, 1 - 2 S. Tongay, 3 H. K. Pal, 3 D. L. Maslov, 3 and A. F. Hebard 3 
1 Institut fur Theorie der kondensierten Materie, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
2 DFG-Center for Functional Nano structures, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
3 Department of Physics, University of Florida, Gainesville, FL 32611, USA 

(Dated: July 2, 2009) 

An interplay between the increase in the number of carriers and the decrease in the scattering 
time is expected to result in a saturation of the in-plane resistivity, p a b, in graphite above room 
temperature. Contrary to this expectation, we observe a pronounced increase in p a b in the interval 
between 300 and 900 K. We provide a theory of this effect based on intervalley scattering of charge 
carriers by high-frequency, graphene-like optical phonons. 

PACS numbers: 81.05.Uw, 72.10-d 



I. INTRODUCTION 

The family of graphite allotropes includes fullerenes, 
carbon nanotubes, graphene, graphite, and diamond. A 
combination of chemical simplicity and diverse physi- 
cal properties, characteristic for these materials, makes 
carbon-based electronics a promising field of research. 
Surprisingly, many properties of the most common mem- 
ber of this group, graphite-such as non-metallic c-axis 
transport and apparent quantization of Hall resistances- 
are yet to be explained. In this work, we focus on in- 
plane transport in graphite at zero magnetic field, which 
is well understood in the degenerate regime, i.e., for 
T < Ep ~ 250 K (we set k B = ft = 1)~ Extending the 
upper limit of the temperature range to 900 K, we ob- 
serve an unexpected pronounced increase in the in-plane 
resistivity. We explain this effect by interaction of carri- 
ers with high-frequency, graphene-like optical phonons. 

The band structure of graphite can be characterized by 
three energy scales X The largest one is set by the hop- 
ping matrix clement between A and B atoms in graphene 
sheets: 70 ~ 3.2 cV. The next two largest ones comprise 
the matrix elements between adjacent Bernal-stacked 
graphene sheets: 71 w 0.35 eV (for the vertical bond) 
and 73 « 0.3 eV (for the slanted bond). Next comes a 
large number of matrix elements between next-to-nearest 
neighbors, which are known to less accuracy, but are gen- 
erally believed not to exceed several tens of meV. The 
major difference between graphite and graphene, i.e., an 
overlap of the conduction and valence band which gives 
rise to a small (»3x 10 18 cm~ 3 at T — > 0) but finite 
carrier concentration, is due to hopping between next-to- 
nearest planes. The band overlap and the Fermi energy 
are on the order of this matrix element. Prior transport 
measurements of graphite were performed at or below 
300 K, where it behaves as a compensated semi-metal 
with finite Fermi energy. There is, however, a very in- 
teresting but hitherto unexplored regime of temperatures 
Ep < T < 71 « 4060 K. This regime can be viewed as a 
critical region of the quantum phase transition between 
a semi-metal with finite band overlap and a semicon- 
ductor with finite band gap. Since next-to-nearest plane 
couplings are irrelevant in this regime, graphite can be 



thought of as a stack of graphene bi-layers.— We dub this 
regime as "bi-layer graphite" (BLGT) . If slanted hopping 
(73) is also neglected, BLGT is described by a simple 
(and historically the first) Wallace model of graphite, 3 
which contains only two hoppings: 70 and 71. The en- 
ergy spectrum in this model consists of two electron and 
two hole branches 



k = ±7ir±J 7 fr 2 + 7 2 |5 k | 



(i) 



where T = cos (k z c/2) , c is the c-axis lattice constant 
and Sk is the structure factor for hopping between in- 
equivalent (A and B) atoms: 



Sk = e ik*a/V3 + 2e- lk ^ a/2 ^cos(ak y /2). 



(2) 



Near the K and K' points of the Brillouin zone, 705k ~ 
vo£k = vq (k x + iky), where vq is the Dirac velocity of a 
single graphene layer. The density of states in BLGT is 
energy-independent v = 167i/w 2 c up to 0(s/ji) terms 



(Ref. [3|). Consequently, the number density of charge 
carriers increases linearly with T and the in-plane con- 
ductivity scales linearly with Tt (Ref. 0) 



lh = (41n2/7r) (e 2 /c) Tt. 



(3) 



At high temperatures, when scattering is predominantly 
due to phonons, one expects 1/t to scale linearly with T 
and, consequently, a a b to be T independent. The con- 
ductivity measured previously up to 300K does indeed 
show a tendency to saturation^ in accord with this ex- 
pectation. However, a different behavior is observed for 
T > 300K. 



II. EXPERIMENT 

We measured the in-plane resistivity of highly oriented 
pyrolytic graphite (HOPG) from 290 K up to 900 K in 
a sealed oven and from 1.7 K up to 310 K in Physical 
Property Measurement System (PPMS) using LR700 16 
Hz AC resistance bridge in four terminal contact con- 
figuration. The samples were cleaved before measure- 
ments. Reproducibility of the results was checked by su- 
perimposing the temperature-dependent resistivities of 
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two separate samples (Fig. [T] inset). All HOPG samples 
were identified to have 0.5 degree mosaic spread deter- 
mined by X-ray rocking curve measurements. Gold con- 
tacts to the sample were made using silver or graphite 
paint. Reproducibility of the measurements was checked 
by sweeping the temperature up and down in the range 
1.7 K < T < 900 K. High temperature measurements 
were performed in a scaled oven with nitrogen flowing 
gas. Temperature was recorded by a type J thermocouple 
located ~ 1mm from the sample to minimize temperature 
lagging effects. Each data point was taken after reaching 
temperature stability. The room temperature resistivity 
was measured to be p^° K ~ 32 ptt-cm before and after 
annealing, consistent with the existing literature values — 
Since p^l 0K remains at the same value before and after 
annealing, we conclude that adsorption/desorption of im- 
purities or doping does not occur up to 900 K. In addi- 
tion, X-ray photoelectron spectroscopy (XPS, Fig. 
and Auger electron spectroscopy (AES, Fig. [2b) were 
performed on different samples and contamination was 
not observed. The XPS and AES spectra displayed re- 
spectively a characteristic Cls peak at 284.6 cV and a 
C auger peak at 271.8 eV. The shape and position of 
the XPS Cls peak was unaltered indicating that C re- 
mained in the same chemical state. When relating the 
measured resistance to resistivity, the lattice constant 
was assumed to be constant. We estimate the change 
of the c-axis lattice constant to be 1% at the highest 
temperature measuredi^ 

The experimental results for p a b{T) are presented in 
Fig. [2 as points. Note that the tendency to saturation 
pronounced at T ~ 300 K is superseded by a rapid in- 
crease which continues unabated up to the highest tem- 
perature measured. The dashed line in the top panel 
shows the theoretical prediction for p a b (T) , calculated 
for 1/r = 1 /to + aT and for a realistic energy spectrum 
of carriers. While this model describes the experiment 
at low temperatures, it fails completely for T > 300 K. 
A slow increase in the theoretical curve, which amounts 
only for a ~ 12% increase of p a b from 300 to 900 K, is 
due to corrections of order T/71 to Eq. ©. 



III. THEORY 



To explain the observed T-dependence of p a b, we first 
notice that the highest resistivity measured, p^ 5K ~ 
62.5 pfl ■ cm, corresponds to 2 M7 per graphene sheet, 
i.e., well below the resistance quantum, h/e 2 « 25 kf2. 
Therefore, the Boltzmann equation should provide an ad- 
equate description of transport in the entire temperature 
interval. To explain the data, we thus need to invoke 
a new scattering mechanism, with 1/r increasing faster 
than T. 
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FIG. 1: (Color online). Measured temperature dependence of 
the in-plane resistivity of HOPG for warming (filled squares) 
and cooling (blank squares) temperature sweeps. Dashed: 
theoretical prediction for p a t(T) in the model containing scat- 
tering at impurities and soft phonons. Solid: fit using the 
model containing scattering at impurities, soft phonons, and 
intervalley scattering at hard in-plane optical phonons [Eq. 
(|17[) ]. Inset: overlap of the data for two different samples. 
The vertical scales were shifted for clarity. 



A. Qualitative picture of transport 

Before proceeding with a discussion of such mecha- 
nisms, it is instructive to develop a more intuitive picture 
of transport in BLGT. To this end, we replace the Wal- 
lace spectrum by a massive (Galilean) form obtained by 

expanding Eq. (fT]) in fen 



Ik 2 + k 2 and keeping only 
the degenerate branches of electrons and holes: 



±k? ] /2m ll (k z ), 



(4) 



where mil (k z ) = Tji/vq is the fc z -dependent mass of the 
in-plane motion. Evaluating the in-plane conductivity as 



Cab (T) 



4e 2 J d 3 k (-dfi/dej wjjr (e k , T) /(2tt) 



(5) 



with spectrum given by Eq. ([?]), r — const, and = 
1/ (exp (ek/T) + 1), we reproduce Eq. ([3]). Typical mo- 
menta contributing to a a b are k z ~ 1/c and hi ~ 
kx = y/2fh\\T, where fh\\ = m\\ (k z = 0) = Ji/vq. Al- 
though expansion in breaks down near the H points 
(k z = ±7r/c), where T vanishes and the spectrum is 
Dirac-like, the contribution of Dirac fermions to <7 a b is 
small in proportion to the volume of the Brillouin zone 
they occupyii Therefore, a typical carrier (in zero mag- 
netic field) in BLGT is massive rather than Dirac-likc 
and the isoenergetic surfaces arc corrugated cylinders 
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FIG. 2: (Color online) . (a) X-ray photoelectron spectroscopy 
(XPS) spectra before and after annealing, (b) Auger electron 
spectroscopy (AES) spectra before and after annealing. In- 
sets: magnified carbon peaks. 



centered near the K points. As in the case of bi-layer 
graphene^ 73 hopping (responsible for trigonal warping 
of the isoenergetic surfaces) leads to a linear- in-fc| | term in 
the energy spectrum, which is smaller than the quadratic 
term for energies > 7173/70 ~ 30 K. Since we are not 
interested here in special effects arising solely from trig- 
onal warping, e.g., longitudinal magnetoresistance, this 
term can be safely neglected. 

Having this simple picture in mind, we now turn to a 
discussion of various scattering mechanisms. 



B. Scattering mechanisms 

1. Electron-hole interaction 

We start with the electron-hole interaction. In contrast 
to the electron-electron interaction, this mechanism gives 
rise to finite resistivity in a compensated semi-metal even 
in the absence of Umklapp processes For T <C Ep, 
we have a usual Fermi- liquid behavior l/r c _h oc T 2 . The 
T 2 -behavior of p a (, is indeed observed in graphite be- 
low 5 K i 11 ! 12 However, BLGT is not a Fermi liquid but 
rather a non-degenerate electron-hole plasma with fixed 
(and equal to zero) chemical potential. To estimate the 
strength of Coulomb interaction, we calculate the r s pa- 
rameter, i.e. the ratio of typical potential and kinetic 



energies 

r s (T) = e 2 /e OQ lE kin , (6) 

where £oo is the background dielectric constant of 
graphite and I = (4tt/3) 1/3 7i" 1/3 oc T~ 1/3 is the typical 
inter-carrier distance. Using an experimentally measured 
values of number density n R* 10 19 cm -3 at T = 300 K 
(Ref.QJ) and 6 W «5 (Ref.QJ) and evaluating the kinetic 
energy as E kin = (7r 2 /121n2) T w 1.2T (Refill!), w e ob- 
tain r s (T = 300 K) w 1.2. As temperature increases, 
r s (T) decreases as T~ 2 / 3 . Therefore, the perturbation 
theory for the electron-hole interaction is reasonably ac- 
curate already at T ~ 300 K and becomes even better 
at higher T. In the Thomas-Fermi model, the screened 

Coulomb potential U (q) = ATre 2 /e 00 (q 2 + q 2 + n 2 ^ is 

isotropic even if the electron spectrum is not; all details of 
the spectrum are encapsulated in the (square of) screen- 
ing wavenumber k 2 = 47re 2 ^/e oc proportional to the den- 
sity of states. In BLGT, k 2 is T-independent and, since 
v is proportional to the (small) in-plane mass, k <C c _1 . 
Also, at not too high temperatures, kr oc y/T <C k. 
Since q\\ cannot exceed the typical electron momentum, 
S| I i5 *C K and U (q) almost docs not depend on q^: 
U (q) w 47re 2 /eoo {q 2 z + n 2 ) . Using this form of U (q) 
and the spectrum from Eq. ((4|), wc obtain from the Fermi 
Golden Rule 

-~^T. (7) 

T c _ h £Kt>o 

This linear scaling persists even at higher temperatures 
(T > 71), in the single-layer limit.— Therefore, electron- 
hole interaction does not provide an explanation of the 
experiment. 

2. Electron-phonon interaction 

We now turn to the electron-phonon scattering. With 
4 atoms per unit cell, graphite has 3 acoustic and 9 op- 
tical phonon modes. The phonon spectrum consists of 
two groups of modes: "hard" and "soft".— The charac- 
teristic energy scale of hard modes, which are present al- 
ready in graphene, is ~ 0.1 eV. Soft modes, with typical 
energies of order 10 meV, arise from weak coupling be- 
tween graphene sheets. For the temperatures of interest 
(T > 300 K), all soft modes are in the classical regime, 
in which the occupation number and, thus, the scatter- 
ing rate scale linearly with T. Although hard acoustic 
modes are still below their Debye temperatures, they 
are also in the classical regime. For example, typical in- 
plane phonon momenta involved in scattering at a hard, 
graphene- like acoustic mode with dispersion uja = Sa6?|| 
are q\\ ~ fen. The corresponding frequencies uja ~ Sa&fc|| 
are smaller than T as long as T > "T,||S 2 fc ~ 1 K. There- 
fore, all soft modes as well as hard acoustic modes lead 
to linear scaling of the scattering rate t~ ^ = aT. 
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The remaining hard modes are graphenc-like optical 
phonons, e.g., the longitudinal optical (i?2g) mode with 
frequency uj ~ 0.17 eV at the T point.— ^ For T < lu , 
scattering at these modes leads to an exponential growth 
of the resistivity with temperature. In the remainder of 
the paper, we will show that this mechanism is capable 
of explaining the experiment. 

Even if only hard phonons are taken into account, a 
real picture of the electron-phonon interaction in graphite 
is rather complicated, as both inter- and intra- valley scat- 
tering on a number of modes are involved. Since our goal 
is just to obtain the temperature dependence of the re- 
sistivity, we will construct a simplified model , 21 > 22 taking 
into account only inter-valley scattering due to one opti- 
cal mode. We start from a tight-binding Hamiltonian 



Hi, H ± 



#graphite — ^jyf jj* J ■ (8) 

Here, #i describes hopping within graphene sheets 



#11 = 



\7o5£ 



7o Sk 

70& 



(9) 



where 5k = icos(-\/3k x a/2)cos(k v a/2) + 2cos(k y a) is 
the structure factor for in-plane hopping between next- 
to-nearest neighbors (AA and BE). As we will show 
shortly, the electron-phonon interaction is dominated 
by the coupling between the diagonal part of #|| and 
phonons. Near the K and K' points, 5 f» —3 + 3fcjja 2 /4. 

The momentum-dependent part of 5 changes the spec- 
trum only at large (~ 1/a) fen and, therefore, can be 
neglected. Hopping in the c-direction is described by 



#i = 



271 cos(fc 2 c/2) 











(10) 



The spectrum of Hamiltonian ||SJ) is given by Eq.([T]), 
while the eigenvectors are represented by a spinor Ak = 
C (±£k/£k, ±1, £k/£k, 1) near the K point and by A£ 
near the K' point. The elements of Ak are the am- 
plitudes of finding a charge carrier on one of the four 
atoms (A, B, A, B) of the graphite unit cell. It is impor- 
tant to notice a difference between Ak and the graphene 
spinors /ik = (l,±£k/£k) ■ Near the K points, £k = |£k| 
in graphene and, therefore, A and B atoms are occupied 
with equal probabilities. If the optical phonon frequency 
is much smaller than 71 , both the initial and final states 
of a scattering process are described by the spectrum in 
Eq. ((4]). For this spectrum, the amplitudes on A and A 

atoms are small as |ek| / |£k| ~ k\\JfhuVQ ~ (T/71) <C 
lf2S Neglecting these amplitudes, we replace the eigen- 
vectors by A k = (0, ±1, 0, 1) / V2. In reality, Wo/71 « 0.5 
and corrections to the results following from this approx- 
imation would be important in a more detailed theory. 

In the simplest model, phonons modulate hopping am- 
plitudes by stretching the corresponding bonds. The 



Hamiltonian of this interaction has the same structure 
as that for an ideal lattice 



#e-ph — 



f! 



(ii) 



where tilde denote corrections to hopping matrix ele- 
ments due to lattice distortions. If #1 does not have diag- 
onal elements, the matrix element of # e _ p h between the 
spinors A k vanishes. Therefore, the electron-phonon in- 
teraction appears only in the diagonal elements, describ- 
ing modulation of hopping between the next-to-nearest 
neighbors. Expanding Hamiltonian (jHJ) in atomic dis- 
placements, one finds 



# 



AA 
e— ph 



sin aq x + 2 sin (aq x /2) cos (^/3aq v /'. 
F y = 2^3 sin (aV3q y /2) cos (aq x /2) , (12) 



where UA,q is the displacement operator for A atoms in 
momentum space. Displacements in the c-direction, be- 
ing perpendicular to the plane, do not change (to linear 
order) the distance between adjacent A atoms, and are 
therefore uncoupled from fermions. Since the formfactors 
F x and F v are small for small q, processes with larger q 
have higher probabilities. For this reason, we will focus 
on inter-valley scattering between K and K' points, cor- 
responding to q = qo = 2n/a (— 1/V3, 1/3) . 

Further expansion in displacements generates higher- 
order vertices of the electron-phonon interaction which 
account for anharmonic effects. Although anharmonic- 
ity is important in a monolayer graphene,— it is weak 
in bulk graphite. For example, the temperature varia- 
tion of the c-axis thermal expansion coefficient amounts 
only to 6% in the temperature interval from 273 to 
1000 K.— Therefore, we neglect the anharmonic effects 
for now but will return to this point when discussing 
the experiment. Employing Eq. (|12p . we find the ma- 

<5k. 



trix element 



.ki+qTV i/2/rxjqL 3 , where 
T q = (c?7o/ da) d q • F is the deformation potential for 
optical phonons, d q is the polarization vector, L is a sys- 
tem size, and p is the mass density. The corresponding 
scattering time can be evaluated as^ 

. (13) 

where the transition rates of emission and absorption of 
phonons, respectively, are 

2^4Vk 2 ! 2 (^+5 ± \) S ki-ei* T <"q) - 

(14) 

and iV q is the Bose function. Neglecting the dispersion of 
the optical mode and using the simplified electron spec- 
trum from Eq. (U), we find for the inter- valley scattering 



5 



rate 

-I, rp\ _ -i coth (lj /2T) cosh 2 ( £k /2T) 

iv 1 k ' J cosh 2 ( ek /2T) + sinh 2 ( Wo /2T)' 1 j 

where all model-dependent details of the electron-phonon 
interaction are incorporated into f _1 . For a thermal elec- 
tron (ejj ~ T) the scattering rate behaves as exp (—loq/T) 
for T <C wo- For T ^> wo, scattering at this mode crosses 
over into the classical regime and Tj" 1 ~ T/luq?. 

At sufficiently high temperatures, where intervalley 
scattering is the dominant mechanism, the conductivity 
is obtained by substituting ^(e^T) into Eq. (|5j) 

a™ = [(4 In 2 - 1) /3tt] (e 2 /c) Tf exp (w /T) , (16) 

C. Comparison to experiment 

Based on the result for the conductivity for intervalley 
scattering [Eq. lfl6|)]. we fit the observed p a b(T) by the 
following formula 

c / 1 m \ 1 C 1 / CJo \ / n 

Pab = - -+aT - + -— r exp (--£ , 17 

where a = 2 (4 In 2 - 1) /3tt r* 0.376 and e* ee 
c / d 3 fcw 2 | (-dfg/dek) /2tt 3 . The first term in Eq. fT7|) 
accounts mostly for the low-T behavior of p a b, when 
scattering at impurities (1/tq) and soft phonons (aT) 
dominate transport. When calculating e* , we accounted 
for hopping between next-to-nearest planes, described by 
72 ~ 0.02 eV. For finite 72, the Fermi energy is finite, 
e* ~ Ep at T — > 0, and the first term in Eq. (fTTj) goes to 
a finite value at T — » 0. The second term is due to inter- 
valley scattering. Equation (fTT)) contains four fitting pa- 
rameters: To, a, f , and ujq. The results of the fit are shown 
in Fig. [TJ The fitting parameters are To = 6.29 x 10~ 12 s, 
a = 0.09, f = 1.4 x 10~ 14 s, and w = 0.22 eV. The values 
of To and a are in reasonable agreement with those found 
previously^ The frequency ujq is somewhat higher but 
still close to the frequency of the £?2g modei 17 ' 18 A rather 
short nominal time f indicates strong coupling between 
electrons and optical phonons in graphite. 

For completeness, we note that for T above 
the Debye frequency, the multi-phonon processes 



modify the scattering rate as 1/t = 1/tq + 
a(T+{T/T 2 ) 2 + {T/T 3 ) 3 + 0{T 4 )), where T 2 , 3 corre- 
spond to an energy scale at which anharmonicity becomes 
strong. The data can be fitted by T 2 « T 3 « 10 3 K, which 
is well below the scale of 10 4 K at which anharmonicity 
becomes strong in other physical properties, such as the 
c-axis thermal expansion coefficient and elastic moduli^ 
This reinforces our conclusion that anharmonicity is not 
important in transport for T < 1000K. 

IV. CONCLUSIONS 

To summarize, we have studied, both experimentally 
and theoretically, the in-plane resistivity of HOPG. We 
found that its temperature dependence is determined by 
a competition between those of the carrier number den- 
sity, n(T), and of the scattering rate 1/t. At tempera- 
tures below 50 K, the number density is practically in- 
dependent of the temperature, while the the scattering 
rate increases with the temperature; as a result, the re- 
sistivity increases with T. At temperatures comparable 
to the Fermi energy, the increase in n(T) almost com- 
pensates for that in 1/t, leading to a quasi-saturation of 
Pab at T ~ 300 K. However, full saturation never occurs 
because, as the temperature increases further, scattering 
off hard optical phonons, characterized by an exponential 
increase of 1/t with T, becomes important. This results 
in a further increase of p a b with T. 
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